Source and mask optimization

ABSTRACT

An illumination source is optimized by changing the intensity and shape of the illumination source to form an image in the image plane that maximizes the minimum ILS at user selected fragmentation points while forcing the intensity at the fragmentation points to be within a small intensity range. An optimum mask may be determined by changing the magnitude and phase of the diffraction orders to form an image in the image plane that maximizes the minimum ILS at user selected fragmentation points while forcing the intensity at the fragmentation points to be within a small intensity range. An optimal mask may be used to create a CPL mask by assigning areas of minimum transmission in an optimum transmission mask a −1, and areas of maximum transmission a +1. Primitive rectangles having a size set to a minimum feature size of a mask maker are assigned to the located minimum and maximum transmission areas and centered at a desired location. The edges of the primitive rectangle are varied to match optimal diffraction orders O(m,n). The optimal CPL mask O CPL (X,Y) is then formed.

[0001] The present invention claims priority from U.S. provisional patent application No. 60/458,365, filed on Mar. 31, 2003, entitled “SOURCE AND MASK OPTIMIZATION.” U.S. provisional patent application No. 60/458,365, filed on Mar. 31, 2003 is herein incorporated by reference in its entirety, including an incorporation of a power point display titled “Source And Mask Optimization For Memory” attached thereto and incorporated by reference in its entirety therein, a power point display titled “Illumination Optimization” also attached thereto and incorporated by reference in its entirety therein, and color drawings labeled FIGS. 1-29 attached thereto.

FIELD OF THE INVENTION

[0002] The field of the invention relates generally to a method and program product for optimizing illumination source and mask features for microlithography.

BACKGROUND

[0003] Lithographic apparatus can be used, for example, in the manufacture of integrated circuits (ICs). In such a case, the photolithographic mask may contain a circuit pattern corresponding to an individual layer of the IC, and this pattern can be imaged onto a target portion (e.g. comprising one or more dies) on a substrate (silicon wafer) that has been coated with a layer of radiation-sensitive material (resist). In general, a single wafer will contain a whole network of adjacent target portions that are successively irradiated via the projection system, one at a time. In one type of lithographic projection apparatus, each target portion is irradiated by exposing the entire mask pattern onto the target portion in one go; such an apparatus is commonly referred to as a wafer stepper. In an alternative apparatus, commonly referred to as a step-and-scan apparatus, each target portion is irradiated by progressively scanning the mask pattern under the projection beam in a given reference direction (the “scanning” direction) while synchronously scanning the substrate table parallel or anti-parallel to this direction. Since, in general, the projection system will have a magnification factor M (generally <1), the speed V at which the substrate table is scanned will be a factor M times that at which the mask table is scanned. More information with regard to lithographic devices as described herein can be gleaned, for example, from U.S. Pat. No. 6,046,792, incorporated herein by reference.

[0004] In a manufacturing process using a lithographic projection apparatus, a mask pattern is imaged onto a substrate that is at least partially covered by a layer of radiation-sensitive material (resist). Prior to this imaging step, the substrate may undergo various procedures, such as priming, resist coating and a soft bake. After exposure, the substrate may be subjected to other procedures, such as a post-exposure bake (PEB), development, a hard bake and measurement/inspection of the imaged features. This array of procedures is used as a basis to pattern an individual layer of a device, e.g. an IC. Such a patterned layer may then undergo various processes such as etching, ion-implantation (doping), metallization, oxidation, chemo-mechanical polishing, etc., all intended to finish off an individual layer. If several layers are required, then the whole procedure, or a variant thereof, will have to be repeated for each new layer. Eventually, an array of devices will be present on the substrate (wafer). These devices are then separated from one another by a technique such as dicing or sawing, whence the individual devices can be mounted on a carrier, connected to pins, etc. Further information regarding such processes can be obtained, for example, from the book Microchip Fabrication: A Practical Guide to Semiconductor Processing, Third Edition, by Peter van Zant, McGraw Hill Publishing Co., 1997, ISBN 0-07-067250-4, incorporated herein by reference.

[0005] For the sake of simplicity, the projection system may hereinafter be referred to as the “lens”; however, this term should be broadly interpreted as encompassing various types of projection systems, including refractive optics, reflective optics, and catadioptric systems, for example. The radiation system may also include components operating according to any of these design types for directing, shaping or controlling the projection beam of radiation, and such components may also be referred to below, collectively or singularly, as a “lens”. Further, the lithographic apparatus may be of a type having two or more substrate tables (and/or two or more mask tables). In such “multiple stage” devices the additional tables may be used in parallel, or preparatory steps may be carried out on one or more tables while one or more other tables are being used for exposures. Twin stage lithographic apparatus are described, for example, in U.S. Pat. No. 5,969,441 and WO 98/40791, incorporated herein by reference.

[0006] The photolithographic masks referred to above comprise geometric patterns corresponding to the circuit components to be integrated onto a silicon wafer. The patterns used to create such masks are generated utilizing CAD (computer-aided design) programs, this process is often being referred to as EDA (electronic design automation). Most CAD programs follow a set of predetermined design rules in order to create functional masks. These rules are set by processing and design limitations. For example, design rules define the space tolerance between circuit devices (such as gates, capacitors, etc.) or interconnect lines, so as to ensure that the circuit devices or lines do not interact with one another in an undesirable way. The design rule limitations are typically referred to as “critical dimensions” (CD). A critical dimension of a circuit can be defined as the smallest width of a line or hole or the smallest space between two lines or two holes. Thus, the CD determines the overall size and density of the designed circuit.

[0007] Various techniques exist to achieve illumination optimization for photolithography. Various mask optimization techniques also have been known. However, currently illumination optimization and mask optimization are not generally linked. U.S. Pat. No. 6,563,566 to Rosenbluth et al. discloses to perform illumination optimization and mask optimization through a series of calculations which attempt to linearize the optimization of the mask transmission. Rosenbluth discloses to maximize the minimum NILS (normalized image log slope) and to select various constraints to be used in the calculations. Rosenbluth also recognizes that the calculations may be limited relying on the symmetry of a mask. However, the linearization of the mask transmission used by Rosenbluth requires using several approximations in the calculations, instead of the actual imaging equations themselves, which produce errors in implementing a mask to form a desired image. The linearization of the mask transmission also requires the use of a significant number of variables, which requires significant computation time to perform the calculations.

[0008] As logic feature sizes decrease, there is a need to provide mask implementations that precisely form a desired image with minimum computational time.

SUMMARY OF THE INVENTION

[0009] In accordance with the present invention, a method for optimizing an illumination source for a mask illumination may comprise the steps of: providing illumination from an illumination source to a plurality of source points and a predetermined mask pattern; selecting fragmentation points in an image plane of an image formed by the illumination provided to the predetermined mask pattern; determining an intensity and image log slope of illumination at each fragmentation point; determining an optimal illumination source as an illumination source which maximizes the image log slope at the selected fragmentation points and has an intensity within a predetermined range.

[0010] In accordance with the present invention, a method for determining an optimal mask may comprise the steps of: determining optimum diffraction orders of an ideal mask; obtaining an optimal transmission mask based on the optimized diffraction orders of the ideal mask; and determining an optimal mask based on the optimal transmission mask, wherein the optimum diffraction orders of the ideal mask are determined by determining a magnitude and phase of diffraction orders which form an image in an image plane which maximizes the minimum illumination log slope at user selected fragmentation points while forcing an intensity of illumination at the fragmentation points to be within a predetermined range.

[0011] In accordance with the present invention, a method of obtaining an optimum source and an optimum mask may comprise the steps of: providing illumination from an illumination source to a plurality of source points and a predetermined mask pattern; selecting fragmentation points in an image plane of an image formed by the illumination provided to the predetermined mask pattern; determining an intensity and image log slope of illumination at each fragmentation point; and simultaneously changing the intensity and shape of the illumination source and the magnitude and phase of diffraction orders of the mask to form an image in the image plane that maximizes the minimum image log slop at the fragmentation points while forcing the intensity at the fragmentation points to be within a predetermined intensity range.

[0012] In accordance with the present invention, a method of optimizing a placement of transmission and phase shifting features on a mask may comprise the steps of: obtaining optimal mask transmission characteristics based on optimum diffraction orders of the mask; locating areas of maximum transmission and minimum transmission; assigning a primitive area as an area centered on an area of maximum transmission or minimum transmission; varying edges of each primitive area to match optimal diffraction orders, wherein each primitive area has a minimum size which is substantially equal to a minimum feature size of the mask.

[0013] In a method of the present invention, the step of obtaining optimal mask transmission characteristics may include a step of determining horizontal diffraction orders of an optimum mask, wherein the number of horizontal diffraction orders is determined according to the equation: $m = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{x}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}$

[0014] where m is the number of horizontal diffraction orders; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.

[0015] In a method of the present invention, the step of obtaining optimal mask transmission characteristics may include a step of determining vertical diffraction orders of an optimum mask, wherein the number of vertical diffraction orders is determined according to the equation: $n = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{y}\left( {\sigma_{\max} + 1} \right)}\quad {NA}}{\lambda} \right\rbrack}} + 1}$

[0016] where n is the number of vertical diffraction orders; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.

[0017] In accordance with the present invention a computer readable medium may contain instructions for a computer to perform a method for optimizing an illumination source for a mask illumination which may comprise the steps: of providing illumination from an illumination source to a plurality of source points and a predetermined mask pattern; selecting fragmentation points in an image plane of an image formed by the illumination provided to the predetermined mask pattern; determining an intensity and image log slope of illumination at each fragmentation point; determining an optimal illumination source as an illumination source which maximizes the image log slope at the selected fragmentation points and has an intensity within a predetermined range.

[0018] In accordance with the present invention a computer readable medium may contain instructions for a computer to cause performance of a method for determining an optimal mask which may comprise the steps of: determining optimum diffraction orders of an ideal mask; obtaining an optimal transmission mask based on the optimized diffraction orders of the ideal mask; and determining an optimal mask based on the optimal transmission mask, wherein the optimum diffraction orders of the ideal mask are determined by determining a magnitude and phase of diffraction orders which form an image in an image plane which maximizes the minimum illumination log slope at user selected fragmentation points while forcing an intensity of illumination at the fragmentation points to be within a predetermined range.

[0019] In accordance with the present invention a computer readable medium may contain instructions for a computer to cause performance of a method of obtaining an optimum source and an optimum mask which may comprise the steps of: providing illumination from an illumination source to a plurality of source points and a predetermined mask pattern; selecting fragmentation points in an image plane of an image formed by the illumination provided to the predetermined mask pattern; determining an intensity and image log slope of illumination at each fragmentation point; and simultaneously changing the intensity and shape of the illumination source and the magnitude and phase of diffraction orders of the mask to form an image in the image plane that maximizes the minimum image log slop at the fragmentation points while forcing the intensity at the fragmentation points to be within a predetermined intensity range.

[0020] In accordance with the present invention a computer readable medium may contain instructions for a computer to cause performance of a method of optimizing a placement of transmission and phase shifting features on a mask which may comprise the steps of: obtaining optimal mask transmission characteristics based on optimum diffraction orders of the mask; locating areas of maximum transmission and minimum transmission; assigning a primitive area as an area centered on an area of maximum transmission or minimum transmission; varying edges of each primitive area to match optimal diffraction orders, wherein each primitive area has a minimum size which is substantially equal to a minimum feature size of the mask.

[0021] In accordance with the present invention a computer readable medium may contain instructions for a computer to cause optimizing a placement of transmission and phase shifting features on a mask comprising the steps of: obtaining optimal mask transmission characteristics; locating areas of minimum transmission; assigning a primitive area as an area centered on an area of minimum transmission; and varying edges of the primitive area to match optimal diffraction orders, wherein the primitive area has a minimum size which is substantially equal to a minimum feature size of the mask.

[0022] A computer readable medium may further contain instructions for a computer to cause the steps of: locating an area of maximum transmission; assigning a transmission primitive area as an area centered on an area of maximum transmission; varying edges of the transmission primitive area to match optimal diffraction orders, wherein the transmission primitive area has a minimum size which is substantially equal to a minimum feature size of the mask.

[0023] In a computer readable medium of the present invention, the step of obtaining optimal mask transmission characteristics may include a step of determining horizontal diffraction orders of an optimum mask, wherein the number of horizontal diffraction orders is determined according to the equation: $m = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{x}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}$

[0024] where m is the number of horizontal diffraction orders; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.

[0025] In a computer readable medium of the present invention, the step of obtaining optimal mask transmission characteristics may include a step of determining vertical diffraction orders of an optimum mask, wherein the number of vertical diffraction orders is determined according to the equation: $n = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{y}\left( {\sigma_{\max} + 1} \right)}\quad {NA}}{\lambda} \right\rbrack}} + 1}$

[0026] where n is the number of vertical diffraction orders; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.

[0027] In accordance with the present invention, an apparatus for optimizing an illumination source for a mask illumination may comprise: an input unit which inputs characteristics of an illumination device; and a processing unit which is configured to change an intensity and shape of an illumination to form an image in an image plane that maximizes the minimum image log slope at user selected fragmentation points.

[0028] In accordance with the present invention, an apparatus for optimizing an mask may comprise: an input unit which inputs a desired image pattern; and a processing unit which is configured to change a magnitude and phase of diffraction orders to form an image in the image plane that maximizes the minimum image log slope at user selected fragmentation points while forcing the intensity at the fragmentation points to be within a predetermined intensity range.

[0029] In accordance with the present invention, an apparatus for obtaining an optimum source and an optimum mask may comprise: an input unit which accepts user inputs; and a processing unit configured to simultaneously change an intensity and shape of an illumination source and change a magnitude and phase of diffraction orders to form an image in an image plane which maximizes a minimum image log slope at user selected fragmentation points while forcing an intensity at the fragmentation points to be within a predetermined intensity range.

[0030] In accordance with the present invention, an apparatus for optimizing a placement of transmission and phase shifting features on a mask comprising: an input unit which inputs characteristics of an illumination device; and a processing unit which is configured to obtain optimal mask transmission characteristics based on optimum diffraction orders of the mask, locate areas of minimum transmission and maximum transmission, assign primitive areas as areas centered on an area of minimum transmission or an area of maximum transmission, and vary edges of the primitive area to match optimal diffraction orders, wherein the primitive areas have a minimum size which is substantially equal to a minimum feature size of the mask.

[0031] In an apparatus of the present invention, the optimal mask transmission characteristics may include horizontal diffraction orders of an optimum mask, and the number of horizontal diffraction orders is determined according to the equation: $m = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{x}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}$

[0032] where m is the number of horizontal diffraction orders; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.

[0033] In an apparatus of the present invention, the optimal mask transmission characteristics may include vertical diffraction orders of an optimum mask, wherein the number of vertical diffraction orders is determined according to the equation: $n = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{y}\left( {\sigma_{\max} + 1} \right)}\quad {NA}}{\lambda} \right\rbrack}} + 1}$

[0034] where n is the number of vertical diffraction orders; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.

[0035] The present invention provides an advantage of a fast computation of mask parameters by limiting the number of mask optimization variables considered. The number of illumination optimization variables may be reduced by relying on the mask symmetry. The number of mask optimization variables may also be reduced by optimizing the diffraction orders of the mask rather than performing computations with the mask transmission. Optimization of the mask diffraction orders is a non-linear process; hence, reducing the number of variables decreases computation time. Furthermore the mask transmission is optimized by performing a non-linear optimization of the diffraction orders followed by a linear optimization of selecting quantized mask transmissions to equal the optimal diffraction orders.

BRIEF DESCRIPTION OF THE DRAWINGS

[0036] The accompanying drawings, which are incorporated in and form a part of the specification, together with the description serve to explain the principles of the invention. In the drawings:

[0037]FIGS. 1A and 1B illustrate an exemplary illumination using Abbe imaging according to the present invention;

[0038]FIG. 2 illustrates an exemplary process for obtaining illumination optimization according to the present invention;

[0039]FIG. 3 provides a graphical representation of the generation of transmission cross coefficients (TCC) according to the present invention;

[0040]FIG. 4 illustrates an exemplary process for performing a mask optimization according to the present invention;

[0041]FIGS. 5A and 5B illustrate additional mask constraints that may be considered in optimizing a mask according to the present invention;

[0042]FIG. 6 illustrates an exemplary process for converting an optimal transmission mask to a CPL mask according to the present invention;

[0043]FIG. 7A illustrates an exemplary DRAM mask pattern and FIG. 7B illustrates an optimum illumination source used with the mask pattern of FIG. 7A;

[0044]FIGS. 8A-8C illustrate diffraction orders of various masks to illustrate the diffraction orders of a CPL mask in accordance with the present invention;

[0045]FIGS. 9A-9C illustrate the application of the optimal illumination to an optimal mask (FIG. 9A);

[0046]FIGS. 10A and 10B illustrate aerial image comparisons between a CPL mask and an 8% AttPSM (phase shifted mask).

[0047]FIGS. 11A and 11B illustrate aerial image comparisons between a CPL mask and an 8% AttPSM using hexapole and annular illumination, respectively;

[0048]FIG. 12A illustrates an exemplary “Short Brickwall” pattern and FIG. 12B illustrates and an illumination source which has been optimized to produce the pattern according to the principles of the present invention;

[0049] FIGS. 13A-C illustrate the diffraction orders of the exemplary mask in FIG. 12A;

[0050] FIGS. 14A-D illustrate the use of primitive edges to create an optimal mask according to the principles of the present invention;

[0051]FIGS. 15A and 15B illustrate aerial image comparisons using the Short Brickwall pattern of FIG. 12A in a CPL mask and an AttPSM mask;

[0052]FIGS. 16A and 16B illustrate the trade-off between obtaining an optimal transmission and using a CPL mask;

[0053]FIGS. 17A and 17B illustrate aerial image comparisons between an 8% AttPSM exposed with a dipole illumination in the top row and an annular illumination in the bottom row;

[0054]FIGS. 18A illustrate the principles of the present invention as applied to a rectangular contact array mask and FIG. 18B illustrates an illumination source which has been optimized to illuminate the patter of FIG. 18A in accordance with the present invention;

[0055] FIGS. 19A-C illustrate the diffraction orders of the exemplary mask in FIG. 18A;

[0056] FIGS. 20A-C illustrate the use of primitive edges to create a quantized CPL mask according to the principles of the present invention;

[0057]FIGS. 21A and 21B illustrate an aerial image comparison using the rectangular contact array mask of FIG. 19A;

[0058]FIGS. 22A and 22B illustrate the trade-off between obtaining an optimal transmission and using a CPL mask;

[0059]FIGS. 23A and 23B illustrate aerial image comparisons between an 8% AttPSM exposed with dipole illumination in the top row and annular illumination in the bottom row;

[0060]FIGS. 24A illustrates a Staggered Rectangular Contact Array and FIG. 24B illustrates and an illumination source which has been optimized to produce the pattern according to the principles of the present invention;

[0061] FIGS. 25A-C illustrate the diffraction orders of the exemplary mask in FIG. 24A;

[0062] FIGS. 26A-C illustrate the use of primitive edges to create a quantized CPL mask according to the principles of the present invention;

[0063]FIGS. 27A and 27B illustrate aerial image comparisons using a CPL mask and an AttPSM mask;

[0064]FIGS. 28A and 28B illustrate the trade-off between obtaining an optimal transmission and using a CPL mask;

[0065]FIGS. 29A and 29B illustrate aerial image comparisons between an 8% AttPSM exposed with quad in the top row and annular illumination in the bottom row;

[0066]FIG. 30 schematically depicts a lithographic projection apparatus suitable for use with a mask designed with the aid of the current invention; and

[0067]FIG. 31 illustrates an exemplary mask optimization processing unit in accordance with the present invention.

DETAILED DESCRIPTION OF THE INVENTION

[0068] Disclosed is a method for optimizing the illumination source and the mask for creating a desired pattern in the image plane. In a preferred implementation, an image is formed in the image plane with the highest image log slope (ILS) in the optimization process at user selected fragmentation points. The image may be optimized by changing the intensity and shape of the illumination source and by changing the magnitude and phase of the diffraction orders. In a preferred implementation, the illumination source is first optimized and then the mask diffraction orders are optimized; however, the illumination source and mask diffraction orders can be optimized simultaneously.

[0069] Since the ILS is a non-linear quantity, the optimization of the illumination source and the mask is a non-linear optimization. Those of skill in the art appreciate that in a non-linear optimization the computational time is a function (e.g. a geometric function) of the number of variables. Therefore to speed the computational time, the number of variables must be minimized. In the illumination source optimization, in accordance with the invention, the number of variables may be reduced by taking advantage of the symmetry of mask. For example, if the mask is symmetric with respect to the vertical and horizontal axes, the illumination source will be symmetric with respect to the vertical and horizontal axes, allowing optimization to be achieved using a quarter of the illumination source.

[0070] In the mask optimization in accordance with the invention, the number of variables may be reduced by performing an optimization of the diffraction orders in the spatial frequency domain. The mask optimization is preferably performed in the spatial frequency domain rather than the spatial domain because the lens (e.g. the entrance pupil of the projection optics) and illumination source limit the number of diffraction orders which are used to form the projected image. In addition to optimizing the ILS at user selected fragmentation points, the shape of the image formed by the mask must match the desired pattern. The matching is preferably done by adding a constraint that intensity at all of the fragmentation points are the same or within a predetermined range of deviation of each other, e.g. a ±2% deviation. Preferably, after optimizing the mask diffraction orders, the optimal mask transmission may be calculated by taking the inverse Fourier transform of the mask diffraction orders.

[0071] In the optimal mask transmission, the electric field transmission has a continuous magnitude from 0 to 1 and a continuous phase from −180 degrees to +180 degrees. Since the mask has a limited number of electric field transmissions, the optimal transmission is preferably quantized by the limited number of transmissions. This quantization is preferably done by selecting quantized mask transmission areas such that the diffraction orders of the quantized mask substantially equal the optimal diffraction orders. Because the Fourier transform is a linear calculation, selecting quantized mask transmission areas such that the quantized mask diffraction orders equal the optimal diffractions is a linear process that can be calculated quickly.

[0072]FIGS. 1A and 1B illustrate an illumination process using Abbe imaging. As illustrated in FIG. 1A, each source point 10 may be illuminated by selectively positioning an illumination source (not shown) at the desired source points 10. The total image intensity is the summation of the intensity from each of the individual source points 10. The illumination patterns are real (in a mathematical sense), hence, the illumination must have even symmetry. Preferably, the source points are chosen to enhance the normalized image log slope (NILS) at fragmentation points on the image plane. A fragmentation point is commonly known to be a point on the image plane which is smaller than λ/2NA.

[0073]FIG. 1B illustrates the use of the illumination source with a typical DRAM mask pattern (referred to as a “Brickwall”). FIG. 7A illustrates a “Long Brickwall” pattern that may be used. In FIG. 1B, portions of the resulting image are shown from illuminating a Brickwall pattern on a mask of a 190 nm pitch with an illumination of λ/800 with a numerical aperture (NA) at the entrance pupil of the projection optics at 0.8. As can be seen in FIG. 1B, the light areas represent an image intensity that enhances the NILS while the dark areas represent an image intensity that degrades the NILS. The illumination source points that provide the best result, e.g. provide the most enhanced NILS, are preferably selected to optimize the shape of the illumination source.

[0074] An exemplary process for performing the illumination source optimization is illustrated in FIG. 2. The illumination source optimization in this process is preferably linear with non-linear constraints. Preferably, the mask transmission magnitude and phase is optimized according to the present invention.

[0075] As illustrated in step S1 of FIG. 2, the user preferably specifies the selected cell (e.g. region on the desired pattern) and the fragmentation points (x,y) (see FIG. 9A) to be evaluated. In the process, illustrated in step S2, a microprocessor preferably calculates the intensity and the NILS at each illuminator point (α,β) and at each fragmentation point (x,y), i.e., the microprocessor calculates I(α,β; x,y) and NILS(α,β; x,y). As also illustrated in step S3 of FIG. 2, the microprocessor utilizes specifications of the illumination system, such as an illumination system from Zeiss (identified as “Zeiss Specs” in FIG. 2), to perform a Gaussian convolution, to determine the minimum pupil fill (e.g., 10%), the minimum ring width (e.g. 0.2) and to force the intensity to a predetermined value. An optimization process is illustrated in step S4, in which illumination points (α,β) that maximize the minimum NILS at each fragmentation point (x,y) are selected. As illustrated in step S5, the intensity I(α,β; x,y) and NILS(α,β; x,y) at each illumination point and each fragmentation point are preferably summed with the selected illumination points that maximize the minimum NILS at each fragmentation point. The optimal illumination source to maximize the NILS at each fragmentation point (x,y) at a desired intensity is determined, as illustrated by step S6.

[0076] Accordingly, in the preferred implementation of the process in FIG. 2, discussed above, the intensity and shape of the illumination may be changed to form an image in the image plane that maximizes the minimum ILS at user selected fragmentation points while forcing the intensity at the fragmentation points to be within a predetermined intensity range.

[0077]FIG. 3 illustrates an exemplary mask optimization process according to the present invention. A technique known as Hopkins Imaging may be used, in which the log slope is attempted to be maximized by changing the pupil to maximize the NILS. Those of skill in the art will appreciate that Abbe Imaging may also be used. Those of skill in the art appreciate that in Abbe Imaging, an image is created for each point and the images are added up and integrated over the source last. Abbe Imaging is generally considered to be spatially incoherent. Those of skill in the art appreciate that in Hopkins Imaging the integration is over the source first and a transfer function is obtained. It may be easier to obtain a mask optimization once from the transmission cross coefficients (TCC) and describe the entire scanner and stepper optics. For mask optimization, the Eigen values decay rapidly to represent the TCC with a few Eigen functions. This speeds computation time.

[0078]FIG. 3 also illustrates the creation of the TCC using Hopkins imaging. The TCC is an autocorrelation of the illumination pupil with the projection pupil. FIG. 3 illustrates the autocorrection of the illumination pupil centered at (0,0) with the projection pupil centered at $\left( {\frac{m\quad \lambda}{P_{x}{NA}},\frac{n\quad \lambda}{P_{y}{NA}}} \right)$

[0079] and with the complex conjugate of the projection pupil centered at $\left( {{- \frac{p\quad \lambda}{P_{x}{NA}}},{- \frac{q\quad \lambda}{P_{y}{NA}}}} \right),$

[0080] where NA represents the numerical aperture of the projection optics, and λ represents the wavelength of the illumination source.

[0081] In performing the Hopkins Imaging, the integration occurs over the source first to form the image transfer cross coefficients (TCC), illustrated as TCC (m,n,p,q) in FIG. 3. The TCC is an autocorrelation of the illumination pupil with the projection pupil and is a four dimensional (4-D) function. The next step is to diagonalize the TCC to reduce the problem to a sum of two dimensional functions. These two dimensional functions are a set of orthogonal eigenfunctions in which each eigenfunction is weighted by the eigenvalue, i.e., eigenfunctions with higher eigenvalues have a larger impact on the image. These eigenfunctions form a set of image kernels which are used in the calculation of the image in the object plane. The diagonalization operation may be performed by any known functions, such as singular value decomposition used in NTI Nanosurfer, or in MG Calibre. Those of skill in the art will also appreciate that a calibrated MT Kernel can also be used.

[0082]FIG. 4 illustrates an exemplary process for performing a mask optimization according to the present invention to obtain an ideal/optimum mask. Ideal mask transmission optimization is non-linear, however the conversion of the ideal mask transmission to a CPL mask implementation is a linear process. In the process illustrated in FIG. 4, the ideal mask is optimized in the frequency domain to speed convergence since the optimization is non-linear. As illustrated in step S21 of FIG. 4, the user selects the cell (e.g. region on the desired pattern) and fragmentation points (x,y) to maximize the NILS and to minimize and maximize intensity. The microprocessor then calculates the TCC (m,n,p,q) (step S22), diagonolizes the TCC into N kernels λφ(m,n) (step S23), and calculates the image intensity to each Kernel, i (step S24). The calculation of the image intensity to each Kernel, i may be performed according to equation 1.

I _(i)(x,y)=|ℑ⁻¹{λ_(i) |O(m,n) ∘φ_(i)(m,n)|}|²  Eq. 1

[0083] In the preferred implementation, the mask transmission range is chosen for CPL mask optimization. The transmission is allowed to be above 1 or below −1 because of the Gibb's phenomenon for image reconstruction after low pass filtering. The mask transmission range can be further modified for an attenuating phase shift mask (PSM). For an attenuating PSM (AttPSM), the mask transmission manufacturable range becomes ℑ⁻¹{O(m,n)}≦1.25 and ℑ⁻¹{O(m,n)}≧−1.25{square root}{square root over (T)} where T is the transmittance of the phase shift mask.

[0084] The microprocessor also receives optimization constraints that may be used to force intensity in the image plane to a predetermined value, minimize intensity below a predetermined value, maximize intensity above a predetermined value, or conform with a mask manufacturability constraint, as illustrated in step S27. The predetermined value is preferably selected as the intensity which provides the highest image log slope (ILS). An exemplary optimization constraint, which may limit the mask transmission to a manufacturable range, is set forth in equation 2.

|ℑ⁻¹ {O(m, n)}|≦1.25  Eq. 2

[0085] The microprocessor preferably performs an optimization of the mask transmissions by changing the diffraction orders of the mask O(m,n) to maximize the NILS at the fragmentation points (x,y), as illustrated in step S25, such as by changing the magnitude and phase of the diffraction orders. The resulting diffraction orders of the optimization of step S25 is summed with the calculated image intensity of each Kernel, i, as illustrated in step S26, and the ideal optimum mask diffraction orders O(m,n) are provided, as illustrated in step S28. The inverse Fourier transform is then performed to convert the calculations to the spatial domain from the frequency domain, as illustrated in step S29, to obtain the optimal transmission mask o(x,y) in the spatial domain, as illustrated in step S30.

[0086] Additional mask constraints may be considered in optimizing a mask as illustrated in FIGS. 5A and 5B. The diffraction order may be evaluated by a real component and an imaginary component and may be represented by equation 3.

O(m,n)=O*(−m,−n)  Eq. 3

[0087] Equation 3 guarantees that the mask is real in a mathematical sense. A real mask has transmission phases of 0° and 180°.

[0088] The number of real diffraction orders, x, may be characterized by equation 4. $\begin{matrix} {x = {{\frac{\left( {m + 1} \right)}{2}\frac{\left( {n + 1} \right)}{2}} + {\frac{\left( {m - 1} \right)}{2}\frac{\left( {n - 1} \right)}{2}}}} & {{Eq}.\quad 4} \end{matrix}$

[0089] The number of imaginary diffraction orders, y, may be characterized by equation 5. $\begin{matrix} {y = {{\frac{\left( {m + 1} \right)}{2}\frac{\left( {n + 1} \right)}{2}} + {\frac{\left( {m - 1} \right)}{2}\frac{\left( {n - 1} \right)}{2}} - 1}} & {{Eq}.\quad 5} \end{matrix}$

[0090] However, a mask must be real, which restricts the number of orders to optimize to x+y, and the entrance pupil of the projection optics limits the number of diffraction orders that may be used by blocking the highest diffraction orders. Hence, the maximum number of horizontal diffraction orders that may be used, m, may be represented by equation 6. $\begin{matrix} {m = {{{floor}\quad\left\lbrack \frac{{P_{x}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack} + 1}} & {{Eq}.\quad 6} \end{matrix}$

[0091] Where m is the number of horizontal diffraction orders, σ_(max) is a radial extent of the distribution of a beam of light from the illumination source, λ is a wavelength of the illumination source, P_(x) is the pitch of the repetitive cell in the x direction, and NA is a numerical aperture of the entrance pupil of the projection optics.

[0092] The number of vertical diffraction orders that may be used, n, may be represented by equation 7. $\begin{matrix} {n = {{{floor}\quad\left\lbrack \frac{{P_{y}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack} + 1}} & {{Eq}.\quad 7} \end{matrix}$

[0093] where n is the number of vertical diffraction orders, σ_(max) is a radial extent of the distribution of a beam of light from the illumination source, λ is a wavelength of the illumination source, P_(y) is the pitch of the repetitive cell in the y direction, and NA is a numerical aperture of the entrance pupil of the projection optics.

[0094] In the preferred implementation, the definition for x and y of the pitch P is the axes of the Cartesian coordinate system in which the x axis is orthogonal to the y axis. However, the coordinate system can be any linear coordinate system in which any two lines, g₁ and g₂, intersecting at an origin, describe the linear coordinate system, i.e. g₁, and g₂, need not be necessarily orthogonal. In such a non-orthogonal coordinate system P_(x) describes the pitch along the axis g₁ and P_(y) describes the pitch along the axis g₂.

[0095] Accordingly, in the preferred implementation of the process illustrated in FIG. 4, as discussed above, the magnitude and phase of the diffraction orders may be changed to form an image in the image plane that maximizes the minimum ILS at user selected fragmentation points while forcing the intensity at the fragmentation points to be within a predetermined intensity range.

[0096] The illumination source optimization, illustrated in FIG. 2, may be simultaneously performed with the mask optimization, illustrated in FIG. 4. Accordingly, the intensity and shape of the illumination and the magnitude and phase of the diffraction orders may be simultaneously changed to form an image in the image plane that maximizes the minimum ILS at user selected fragmentation points while forcing the intensity at the fragmentation points to be within a predetermined intensity range.

[0097] The ideal optimal transmission mask determined in the exemplary process of FIG. 4 maybe converted to be implemented in an actual mask, such as a CPL mask, as illustrated by the process illustrated in FIG. 6. As shown in step S31 in FIG. 6, an ideal optimal transmission mask, which may be determined according to the process illustrated in FIG. 4, is provided to be converted. In the process, a dark field mask is preferably used to start, as illustrated in step S32. Then areas of minimum transmission are located and assigned a −1 (step S33), and areas of maximum transmission are also located and assigned a +1 (step S34). Primitive rectangles preferably having a size set to a minimum feature size of a mask maker are assigned to the located minimum transmission areas and centered at a desired location (step S35). Likewise, primitive rectangles are assigned for the located areas of maximum transmission and centered (step S36). With the assigned −1 and +1 values, the diffraction orders needed for optimization start at ${- 1}*{{floor}\quad\left\lbrack \frac{{P_{x}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}$

[0098] and end at ${+ 1}*{{{floor}\quad\left\lbrack \frac{{P_{x}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}.}$

[0099] Hence, the maximum number of horizontal diffraction orders that may be used, m, may be represented by the equation: $\begin{matrix} {m = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{x}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}} & {{Eq}.\quad 8} \end{matrix}$

[0100] and the equation: $\begin{matrix} {n = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{y}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}} & {{Eq}.\quad 9} \end{matrix}$

[0101] may represent the maximum number of vertical diffraction orders, n that may be used. As illustrated in step S37, the −1 and +1 rectangle edges may be varied to match optimal diffraction orders O(m,n). The optimal CPL mask O_(CPL)(x,y) may be formed as illustrated in step S38.

[0102] In the process illustrated in FIG. 6, the mask conversion is linear. However, CPL edge effects can also be taken into account by modifying the diffraction order through a perturbational model. In the perturbational model, edges of the mask are replaced with local areas that have transmission of neither 0, +1, or −1. These areas allow a scalar mask to emulate the vector edge effects from a CPL mask. Those of skill in the art will appreciate that many perturbational models may be used, such as those disclosed by J. Tirapu-Azpiroz, E. Yablonovitch, “Fast evaluation of Photomask, Near-Fields in Sub-Wavelength 193 nm Lithography,” Proc of the SPIE, vol. 5377 (2004), and K. Adam, A. Neureuther, “Simplified Models For Edge Transitions In Rigorous Mask Modeling,” Proc. of the SPIE, vol. 4346 (2001), pp. 331-344.

[0103]FIG. 7A illustrates an exemplary DRAM mask pattern, commonly referred to as a Long Brickwall. FIG. 7B illustrates an optimum illumination source for the mask of FIG. 7A obtained according to the principles of the present invention. The illumination source and CPL mask are optimized for λ/800, NA=0.8 and 190 nm pitch (k1=0.306). The source illumination was optimized for the long brickwall using the algorithm outlined in FIG. 2. The illumination poles on the y axis help improve the NILS at the end of the long brick while the illumination poles on the x axis help improve the NILS in between the bricks.

[0104]FIGS. 8A-8C illustrate diffraction orders of various masks. FIG. 8A illustrates the original diffraction orders of the mask utilizing AttPSM. FIG. 8B illustrates the ideal optimal diffraction orders determined according to the process illustrated in FIG. 4 of the present invention. FIG. 8C illustrates the implementation of the optimal diffraction orders with a CPL mask, according to the process illustrated in FIG. 6 of the present invention. The optimal diffraction orders force more energy into the higher orders (±2,0) and (±1, ±3). Those of skill in the art will appreciate that a CPL mask can represent the optimal diffraction orders almost identically.

[0105]FIGS. 9A-9C illustrate the application of the optimal illumination to an optimal mask (FIG. 9A). In FIG. 9A the “*” points identify the NILS optimization points, the “+” points identify the points at which intensity should be maximized, and the “−” points identify the points at which intensity should be minimized. The NILS optimization attempts to make the slope of the intensity at the image plane as large as possible to obtain a high contrast in printing the mask features. The NILS is preferably within a set value which is suitable for the determined NILS optimization points along a mask feature. Those of skill in the art will appreciate that a CPL mask has three phase values, +180, 0 and −180 degrees. The CPL mask illustrated in FIG. 8C is formed by quantizing the optimal mask transmission, to form a quantized CPL mask transmission.

[0106]FIGS. 10A and 10B and 11A and 11B illustrate aerial image comparisons between a CPL mask and an 8% AttPSM (attenuated phase shifted mask). In FIG. 10A the CPL mask was exposed with hexapole illumination. In FIG. 10B the PSM mask was exposed with hexapole illumination. As illustrated in FIGS. 10A and 10B, the contrast and the NILS is much better in between the lines with the CPL mask than the PSM mask. It can be seen that CPL has a production worthy process margin while the 8% solution is less favorable. However, CPL has a little necking which may be correctable by adding more optimization points.

[0107] Also, as illustrated in FIG. 11A, the hexapole illuminator for 8% Att PSM also shows NILS improvement over annular illumination shown in FIG. 11B. However, a larger process margin may be gained by using CPL with the hexapole illuminator.

[0108] Optimizing the mask in the frequency domain limits the size of the optimization problem and speeds convergence. Reconstruction of the optimal CPL mask from the optimal diffraction orders is a linear problem in the frequency domain. Process window improvement may be optimal with the optimal diffraction orders for all structures, which may minimize the effects of focus and exposure variations in the scanner and in the wafer topography. Process window as used herein is the amount of exposure latitude (EL) as a function of depth of focus (DOF). Process window improvement with the optimal CPL may also be optimal. CPL may be used to improve the process window at a low k1 factor. The k1 factor as used herein may be defined by CD*NA/λ, where CD is the critical dimension of the feature to be printed, and λ is the wavelength of the illumination source. Additionally, illumination optimization does not need to run first. Using the present invention the mask can be optimized to an existing OAI (quasar, cquad (a quadrupole illuminator with poles on the Cartesian x and y axes, such as a quasar illuminator rotated 45 degrees), annular illumination). Also, in the present invention, the mask can be optimized for a single or double exposure (1 or 2 masks) or for a 2 tone mask (binary or AttPSM). Most ideal transmission masks can be represented with CPL in accordance with the present invention.

[0109]FIGS. 12A and 12B illustrate an exemplary “Short Brickwall” pattern in which the source and the CPL mask have been optimized according to the principles of the present invention. FIG. 12B illustrates an optimum illumination source for the pattern of FIG. 12A which has been obtained according to the principles of the present invention as illustrated in FIG. 2. The source and the CPL mask were optimized for λ/800 NA=0.8 and 190 nm pitch. Once again, the contrast and the NILS can be observed.

[0110] FIGS. 13A-C illustrate the diffraction orders of the exemplary mask in FIG. 12A. FIG. 13A illustrates the original diffraction order with the AttPSM mask, FIG. 13B illustrates the ideal optimal diffraction orders determined according to the process illustrated in FIG. 4 of the present invention. FIG. 13C illustrates an implementation of the optimal diffraction orders with a CPL mask according to the process illustrated in FIG. 4 of the present invention. Those of skill in the art will appreciate that the optimal diffraction orders force more energy into (±1, ±1) areas. As illustrated by FIGS. 13B and 13C, the CPL mask can represent the optimal diffraction orders almost identically.

[0111] FIGS. 14A-D illustrate the use of primitive edges to create an optimal mask according to the principles of the present invention. FIG. 14A illustrates the exemplary Short Brickwall mask of FIG. 12A. FIG. 14B illustrates the optimal mask transmission of the mask in FIG. 14A. FIG. 14C illustrates using an arcuate modeling technique to map the primitive edges more closely with the optimal mask transmission. FIG. 14D illustrates using primitive rectangles to map the optimal mask transmission. As illustrated by FIGS. 14C and 14D, two CPL implementations may lead to substantially the same diffraction order spectrum. Primitive edges are not critical but using primitives may lead to a smaller figure count and possibly allow an easier inspection of the mask.

[0112]FIGS. 15A-17B illustrate aerial image comparisons using the Short Brickwall of FIG. 12A. In FIG. 15, a CPL mask and an 8% AttPSM mask are both exposed with dipole illumination. As seen in FIGS. 15A and 15B, the contrast and NILS is much better in between the lines and the end of the line with CPL than the PSM mask. The CPL mask has production worthy process margins while the 8% solution does not appear to be as favorable. As also illustrated in FIG. 15A and 15B, the CPL mask maintains the area better than 8% AttPSM.

[0113]FIGS. 16A and 16B illustrates the trade-off between obtaining an optimal transmission and using a CPL mask. As illustrated by FIGS. 16A and 16B, there is almost no difference between optimal transmission and CPL representation. Hence, the CPL mask offers a favorable solution for mask optimization according to the principles of the present invention.

[0114]FIGS. 17A illustrates aerial image comparisons between an 8% AttPSM exposed with dipole illumination and FIG. 17B illustrates using annular illumination. As shown in FIG. 17A and 17B, the NILS is better with dipole illumination in comparison to annular illumination. Annular illumination has better NILS at the end of line (EOL). However, the NILS probably is not large enough with annular illumination to print bricks without bridging.

[0115]FIGS. 18A-23B illustrate the principles of the present invention as applied to a rectangular contact array mask, depicted in FIG. 18A. In FIGS. 18A and 18B, with the source and the CPL Mask optimized for λ/800 NA=0.8 and 190 nm pitch. Once again, the contrast and the NILS can be observed.

[0116] FIGS. 19A-C illustrate the diffraction orders of the exemplary mask in FIG. 18A. FIG. 19A illustrates the original diffraction order with the AttPSM mask, FIG. 19B illustrates the optimal diffraction order determined by the principles of the present invention, and FIG. 19C illustrates the optimal diffraction order implemented with a CPL mask. Those of skill in the art will appreciate that the optimal diffraction orders force more energy into the higher orders (±1, 0), (0, ±), and (0, ±2). As illustrated by FIGS. 19B and 19C, the CPL mask can represent the optimal diffraction orders almost identically.

[0117] FIGS. 20A-C illustrate the use of primitive edges to create a quantized CPL mask according to the principles of the present invention. FIG. 20A illustrates the exemplary Rectangular Contact Array mask of FIG. 19A. FIG. 20B illustrates the optimal mask transmission of the mask in FIG. 20A. FIG. 20C illustrates using primitive rectangles to map the optimal mask transmission to create a quantized CPL mask.

[0118]FIGS. 21A-23B illustrate aerial image comparisons using the Rectangular Contact Array mask of FIG. 19A. In FIG. 21A, a CPL mask and in FIG. 21B an 8% AttPSM mask are both exposed with hexapole illumination. As seen in FIGS. 21A and 21B, the peak intensity and NILS are better with the CPL mask. The CPL mask has production worthy DOF (depth of focus) while the 8% AttPSM solution does not appear to have sufficient DOF to be as favorable. As also illustrated in FIGS. 21A and 21B, the CPL mask maintains the area better than the 8% AttPSM.

[0119]FIGS. 22A and 22B illustrates a trade-off between obtaining an optimal transmission and using a CPL mask. As illustrated by FIGS. 22A and 22B, there appears to be a slightly better NILS with optimal transmission reticle in comparison to CPL.

[0120]FIGS. 23A and 23B illustrates aerial image comparisons between an 8% AttPSM exposed with dipole illumination in the top row and annular illumination in the bottom row. As shown in FIGS. 23A and 23B, the peak intensity and NILS are better with hexapole illumination in comparison to annular illumination. The peak intensity with annular illumination may not be sufficient to print through focus.

[0121]FIGS. 24A-29B illustrate the principles of the present invention as applied to a Staggered Rectangular Contact Array, depicted in FIG. 24A. FIG. 24A and 24B illustrate an exemplary Staggered Rectangular Contact Array mask in which the source and CPL mask have been optimized according to the principles of the present invention. In FIGS. 24A and 24B, the source and CPL Mask were optimized for λ/800, NA=0.8 and 190 nm pitch. Once again, the contrast and the NILS can be observed.

[0122] FIGS. 25A-C illustrate the diffraction orders of the exemplary mask in FIG. 24A. FIG. 25A illustrates the original diffraction order with the AttPSM mask, FIG. 25B illustrates the optimal diffraction order determined by the principles of the present invention, and FIG. 25C illustrates the optimal diffraction order with a CPL implementation. Those of skill in the art will appreciate that the optimal diffraction orders force more energy into (0, +2) and (±1, ±1) areas. As illustrated by FIGS. 25B and 25C, the CPL mask can represent the optimal diffraction orders almost identically.

[0123] FIGS. 26A-C illustrate the use of primitive edges to create a quantized CPL mask according to the principles of the present invention. FIG. 26A illustrates the exemplary Staggered Rectangular Contact Array mask of FIG. 24A. FIG. 26B illustrates the optimal mask transmission of the mask in FIG. 24A. FIG. 26C illustrates using primitive rectangles to map the optimal mask transmission to create a quantized CPL mask. In FIG. 26A-C, 180 degree outriggers may be used at the sides of the contact array.

[0124]FIGS. 27A-29B illustrate aerial image comparisons using the Rectangular Contact Array mask of FIG. 24A. In FIGS. 27A and 27B, a CPL mask and an 8% AttPSM, respectively, mask are both exposed with quad illumination. As seen in FIG. 27A and 27B, the peak intensity and NILS are better with the CPL mask than the PSM mask. The CPL mask also has a greater exposure latitude and DOF in comparison to the 8% AttPSM mask.

[0125]FIGS. 28A and 28B illustrates the trade-off between obtaining an optimal transmission and using a CPL mask. As illustrated by FIGS. 28A and 28B, there is almost no difference between optimal transmission and the CPL mask representation. The CPL mask offers a favorable solution for mask optimization according to the principles of the present invention.

[0126]FIGS. 29A and 29B illustrates aerial image comparisons between an 8% AttPSM exposed with quad illumination in the top row and annular illumination in the bottom row. As shown in FIGS. 29A and 29B, the peak intensity and NILS are better with quad illumination than annular illumination. The peak intensity with annular illumination may not be sufficient to print through focus.

[0127]FIG. 30 schematically depicts a lithographic projection apparatus suitable for use with a mask designed with the aid of the current invention. The apparatus comprises:

[0128] a radiation system Ex, IL, for supplying a projection beam PB of radiation. In this particular case, the radiation system also comprises a radiation source LA;

[0129] a first object table (mask table) MT provided with a mask holder for holding a mask MA (e.g. a reticle), and connected to first positioning means for accurately positioning the mask with respect to item PL;

[0130] a second object table (substrate table) WT provided with a substrate holder for holding a substrate W (e.g. a resist-coated silicon wafer), and connected to second positioning means for accurately positioning the substrate with respect to item PL;

[0131] a projection system (“lens”) PL (e.g. a refractive, catoptric or catadioptric optical system) for imaging an irradiated portion of the mask MA onto a target portion C (e.g. comprising one or more dies) of the substrate W.

[0132] As depicted herein, the apparatus is of a transmissive type (i.e. has a transmissive mask). However, in general, it may also be of a reflective type, for example (with a reflective mask). Alternatively, the apparatus may employ another kind of patterning means as an alternative to the use of a mask; examples include a programmable mirror array or LCD matrix.

[0133] The source LA (e.g. a mercury lamp or excimer laser) produces a beam of radiation. This beam is fed into an illumination system (illuminator) IL, either directly or after having traversed conditioning means, such as a beam expander Ex, for example. The illuminator IL may comprise adjusting means AM for setting the outer and/or inner radial extent (commonly referred to as σ-outer and σ-inner, respectively) of the intensity distribution in the beam. In addition, it will generally comprise various other components, such as an integrator IN and a condenser CO. In this way, the beam PB impinging on the mask MA has a desired uniformity and intensity distribution in its cross-section.

[0134] It should be noted with regard to FIG. 30 that the source LA may be within the housing of the lithographic projection apparatus (as is often the case when the source LA is a mercury lamp, for example), but that it may also be remote from the lithographic projection apparatus, the radiation beam that it produces being led into the apparatus (e.g. with the aid of suitable directing mirrors); this latter scenario is often the case when the source LA is an excimer laser (e.g. based on KrF, ArF or F₂ lasing). The illumination source intensity may also be made with a mirror array or an LCD. The current invention encompasses at least both of these scenarios.

[0135] The beam PB subsequently intercepts the mask MA, which is held on a mask table MT. Having traversed the mask MA, the beam PB passes through the lens PL, which focuses the beam PB onto a target portion C of the substrate W. With the aid of the second positioning means (and interferometric measuring means IF), the substrate table WT can be moved accurately, e.g. so as to position different target portions C in the path of the beam PB. Similarly, the first positioning means can be used to accurately position the mask MA with respect to the path of the beam PB, e.g. after mechanical retrieval of the mask MA from a mask library, or during a scan. In general, movement of the object tables MT, WT will be realized with the aid of a long-stroke module (coarse positioning) and a short-stroke module (fine positioning), which are not explicitly depicted in FIG. 27. However, in the case of a wafer stepper (as opposed to a step-and-scan tool) the mask table MT may just be connected to a short stroke actuator, or may be fixed.

[0136] The depicted tool can be used in two different modes:

[0137] In step mode, the mask table MT is kept essentially stationary, and an entire mask image is projected in one go (i.e. a single “flash”) onto a target portion C. The substrate table WT is then shifted in the x and/or y directions so that a different target portion C can be irradiated by the beam PB;

[0138] In scan mode, essentially the same scenario applies, except that a given target portion C is not exposed in a single “flash”. Instead, the mask table MT is movable in a given direction (the so-called “scan direction”, e.g. the y direction) with a speed ν, so that the projection beam PB is caused to scan over a mask image; concurrently, the substrate table WT is simultaneously moved in the same or opposite direction at a speed V=Mν, in which M is the magnification of the lens PL (typically, M=¼ or ⅕). In this manner, a relatively large target portion C can be exposed, without having to compromise on resolution.

[0139] The concepts disclosed herein may simulate or mathematically model any generic imaging system for imaging sub wavelength features, and may be especially useful with emerging imaging technologies capable of producing wavelengths of an increasingly smaller size. Emerging technologies already in use include EUV (extreme ultra violet) lithography that is capable of producing a 193 nm wavelength with the use of a ArF laser, and even a 157 nm wavelength with the use of a Fluorine laser. Moreover, EUV lithography is capable of producing wavelengths within a range of 20-5 nm by using a synchrotron or by hitting a material (either solid or a plasma) with high energy electrons in order to produce photons within this range. Because most materials are absorptive within this range, illumination may be produced by reflective mirrors with a multi-stack of Molybdenum and Silicon. The multi-stack mirror has a 40 layer pairs of Molybdenum and Silicon where the thickness of each layer is a quarter wavelength. Even smaller wavelengths may be produced with X-ray lithography. Typically, a synchrotron is used to produce an X-ray wavelength. Since most material is absorptive at x-ray wavelengths, a thin piece of absorbing material defines where features would print (positive resist) or not print (negative resist).

[0140] While the concepts disclosed herein may be used for imaging on a substrate such as a silicon wafer, it shall be understood that the disclosed concepts may be used with any type of lithographic imaging systems, e.g., those used for imaging on substrates other than silicon wafers.

[0141] Software functionalities of a computer system involve programming, including executable code, may be used to implement the above described imaging model. The software code is executable by the general-purpose computer. In operation, the code and possibly the associated data records are stored within a general-purpose computer platform. At other times, however, the software may be stored at other locations and/or transported for loading into the appropriate general-purpose computer systems. Hence, the embodiments discussed above involve one or more software products in the form of one or more modules of code carried by at least one machine-readable medium. Execution of such code by a processor of the computer system enables the platform to implement the catalog and/or software downloading functions, in essentially the manner performed in the embodiments discussed and illustrated herein.

[0142] As used herein, terms such as computer or machine “readable medium” refer to any medium that participates in providing instructions to a processor for execution. Such a medium may take many forms, including but not limited to, non-volatile media, volatile media, and transmission media. Non-volatile media include, for example, optical or magnetic disks, such as any of the storage devices in any computer(s) operating as one of the server platform, discussed above. Volatile media include dynamic memory, such as main memory of such a computer platform. Physical transmission media include coaxial cables; copper wire and fiber optics, including the wires that comprise a bus within a computer system. Carrier-wave transmission media can take the form of electric or electromagnetic signals, or acoustic or light waves such as those generated during radio frequency (RF) and infrared (IR) data communications. Common forms of computer-readable media therefore include, for example: a floppy disk, a flexible disk, hard disk, magnetic tape, any other magnetic medium, a CD-ROM, DVD, any other optical medium, less commonly used media such as punch cards, paper tape, any other physical medium with patterns of holes, a RAM, a PROM, and EPROM, a FLASH-EPROM, any other memory chip or cartridge, a carrier wave transporting data or instructions, cables or links transporting such a carrier wave, or any other medium from which a computer can read programming code and/or data. Many of these forms of computer readable media may be involved in carrying one or more sequences of one or more instructions to a processor for execution.

[0143] As illustrated in FIG. 31, an exemplary mask optimization unit may contain a processor 1000 which receives input from an input unit 1003. Processor 1000 may be a conventional microprocessor or may be a specially designed processing unit, such as an EEPROM or EPROM or a fabricated integrated circuit. Input 1003 may be any type of electronic input device, such as a keyboard or a mouse, or may be a memory or internet connection. Processor 1000 preferably retrieves stored protocols from ROM 1002 and RAM 1001, such as protocols to implement the processing illustrated in FIGS. 2-6, and stores information on RAM 1001. The calculated results of processor 1000 may be displayed on display 1004 and may be provided to a mask fabrication unit.

[0144] The invention may be embodied in other specific forms without departing from the spirit or essential characteristics thereof. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. 

what is claimed is:
 1. A method for optimizing an illumination source for a mask illumination comprising the steps of: providing illumination from an illumination source to a plurality of source points and a predetermined mask pattern; selecting fragmentation points in an image plane of an image formed by the illumination provided to the predetermined mask pattern; determining an intensity and image log slope of illumination at each fragmentation point; determining an optimal illumination source as an illumination source which maximizes the image log slope at the selected fragmentation points and has an intensity within a predetermined range.
 2. The method of claim 1, further providing the step of providing optimization constraints on the optimal illumination source.
 3. The method of claim 1, wherein the step of determining the intensity and shape of the illumination source, determines an optimal illumination source as one which forces illumination intensity at the image plane to a predetermined value.
 4. The method of claim 1, wherein the step of determining an optimal illumination source determines an optimal shape of the illumination source.
 5. The method of claim 1, wherein the step of determining an optimal illumination source determines an optimal intensity of the illumination source.
 6. A method for determining an optimal mask comprising the steps of: determining optimum diffraction orders of an ideal mask; obtaining an optimal transmission mask based on the optimized diffraction orders of the ideal mask; and determining an optimal mask based on the optimal transmission mask, wherein the optimum diffraction orders of the ideal mask are determined by determining a magnitude and phase of diffraction orders which form an image in an image plane which maximizes the minimum illumination log slope at user selected fragmentation points while forcing an intensity of illumination at the fragmentation points to be within a predetermined range.
 7. The method of claim 6, wherein the step of obtaining optimal mask transmission characteristics includes a step of determining horizontal diffraction orders of an optimum mask, wherein the number of horizontal diffraction orders is determined according to the equation: $m = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{x}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}$

where m is the number of horizontal diffraction orders; P_(x) is the pitch of the repetitive cell in the x direction; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.
 8. The method of claim 6, wherein the step of obtaining optimal mask transmission characteristics includes a step of determining vertical diffraction orders of an optimum mask, wherein the number of vertical diffraction orders is determined according to the equation $n = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{y}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}$

where n is the number of vertical diffraction orders; P_(y) is the pitch of the repetitive cell in the y direction; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.
 9. The method of claim 6, wherein the step of determining optimum diffraction orders determines optimum diffraction orders in the special frequency domain.
 10. The method of claim 6, wherein the step of determining an optimum mask comprises the steps of: locating areas of maximum transmission and minimum transmission; assigning a primitive area as an area centered on an area of maximum transmission or minimum transmission; varying edges of each primitive area to match optimal diffraction orders, wherein each primitive areas has a minimum size which is substantially equal to a minimum feature size of the mask.
 11. A method of obtaining an optimum source and an optimum mask comprising the steps of: providing illumination from an illumination source to a plurality of source points and a predetermined mask pattern; selecting fragmentation points in an image plane of an image formed by the illumination provided to the predetermined mask pattern; determining an intensity and image log slope of illumination at each fragmentation point; and simultaneously changing the intensity and shape of the illumination source and the magnitude and phase of diffraction orders of the mask to form an image in the image plane that maximizes the minimum image log slop at the fragmentation points while forcing the intensity at the fragmentation points to be within a predetermined intensity range.
 12. A method of optimizing a placement of transmission and phase shifting features on a mask comprising the steps of: obtaining optimal mask transmission characteristics based on optimum diffraction orders of the mask; locating areas of maximum transmission and minimum transmission; assigning a primitive area as an area centered on an area of maximum transmission or minimum transmission; varying edges of each primitive area to match optimal diffraction orders, wherein each primitive area has a minimum size which is substantially equal to a minimum feature size of the mask.
 13. The method of claim 12, wherein the step of obtaining optimal mask transmission characteristics includes a step of determining horizontal diffraction orders of an optimum mask, wherein the number of horizontal diffraction orders is determined according to the equation: $m = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{x}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}$

where m is the number of horizontal diffraction orders; P_(x) is the pitch of the repetitive cell in the x direction; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.
 14. The method of claim 12, wherein the step of obtaining optimal mask transmission characteristics includes a step of determining vertical diffraction orders of an optimum mask, wherein the number of vertical diffraction orders is determined according to the equation $n = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{y}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}$

where n is the number of vertical diffraction orders; P_(y) is the pitch of the repetitive cell in the y direction; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.
 15. The method of claim 12, wherein the mask is a CPL mask.
 16. A computer readable medium containing instructions for a computer to cause optimizing an illumination source for a mask illumination comprising the steps of: providing illumination from an illumination source to a plurality of source points and a predetermined mask pattern; selecting fragmentation points in an image plane of an image formed by the illumination provided to the predetermined mask pattern; determining an intensity and image log slope of illumination at each fragmentation point; determining an optimal illumination source as an illumination source which maximizes the image log slope at the selected fragmentation points and has an intensity within a predetermined range.
 17. The computer readable medium of claim 16, further containing instructions for a computer to cause the step of providing optimization constraints on the optimal illumination source.
 18. The computer readable medium of claim 16, wherein the step of determining the intensity and shape of the illumination source, determines an optimal illumination source as one which forces illumination intensity at the image plane to a predetermined value.
 19. The computer readable medium of claim 16, wherein the step of determining an optimal illumination source determines an optimal shape of the illumination source.
 20. The computer readable medium of claim 16, wherein the step of determining an optimal illumination source determines an optimal intensity of the illumination source.
 21. A computer readable medium containing instructions for a computer to perform a method for determining an optimal mask comprising the steps of: determining optimum diffraction orders of an ideal mask; obtaining an optimal transmission mask based on the optimized diffraction orders of the ideal mask; and determining an optimal mask based on the optimal transmission mask, wherein the optimum diffraction orders of the ideal mask are determined by determining a magnitude and phase of diffraction orders which form an image in an image plane which maximizes the minimum illumination log slope at user selected fragmentation points while forcing an intensity of illumination at the fragmentation points to be within a predetermined range.
 22. The computer readable medium of claim 21, wherein the step of obtaining optimal mask transmission characteristics includes a step of determining horizontal diffraction orders of an optimum mask, wherein the number of horizontal diffraction orders is determined according to the equation: $m = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{x}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}$

where m is the number of horizontal diffraction orders; P_(x) is the pitch of the repetitive cell in the x direction; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.
 23. The computer readable medium of claim 21, wherein the step of obtaining optimal mask transmission characteristics includes a step of determining vertical diffraction orders of an optimum mask, wherein the number of vertical diffraction orders is determined according to the equation $n = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{y}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}$

where n is the number of vertical diffraction orders; P_(y) is the pitch of the repetitive cell in the y direction; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.
 24. The method of claim 21, wherein the step of determining optimum diffraction orders determines optimum diffraction orders in the special frequency domain.
 25. The computer readable medium of claim 21, wherein the step of determining an optimum mask comprises the steps of: locating areas of maximum transmission and minimum transmission; assigning a primitive area as an area centered on an area of maximum transmission or minimum transmission; varying edges of each primitive area to match optimal diffraction orders, wherein each primitive areas has a minimum size which is substantially equal to a minimum feature size of the mask.
 26. A computer readable medium containing instructions for a computer to perform a method of obtaining an optimum source and an optimum mask comprising the steps of: providing illumination from an illumination source to a plurality of source points and a predetermined mask pattern; selecting fragmentation points in an image plane of an image formed by the illumination provided to the predetermined mask pattern; determining an intensity and image log slope of illumination at each fragmentation point; and simultaneously changing the intensity and shape of the illumination source and the magnitude and phase of diffraction orders of the mask to form an image in the image plane that maximizes the minimum image log slop at the fragmentation points while forcing the intensity at the fragmentation points to be within a predetermined intensity range.
 27. A computer readable medium containing instructions for a computer to cause optimizing a placement of transmission and phase shifting features on a mask comprising the steps of: obtaining optimal mask transmission characteristics based on optimum diffraction order of the mask; locating areas of maximum transmission and minimum transmission; assigning a primitive area as an area centered on an area of maximum transmission or minimum transmission; varying edges of each primitive area to match optimal diffraction orders, wherein each primitive area has a minimum size which is substantially equal to a minimum feature size of the mask.
 28. The computer readable medium of claim 27, wherein the mask is a CPL mask.
 29. The computer readable medium of claim 27, wherein the step of obtaining optimal mask transmission characteristics includes a step of determining horizontal diffraction orders of an optimum mask, wherein the number of horizontal diffraction orders is determined according to the equation $m = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{x}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}$

where m is the number of horizontal diffraction orders; P_(x) is the pitch of the repetitive cell in the x direction; λ is a_(wavelength) of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.
 30. The computer readable medium of claim 27, wherein the step of obtaining optimal mask transmission characteristics includes a step of determining vertical diffraction orders of an optimum mask, wherein the number of vertical diffraction orders is determined according to the equation: $n = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{y}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}$

where n is the number of vertical diffraction orders; P_(y) is the pitch of the repetitive cell in the y direction; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.
 31. An apparatus for optimizing an illumination source for a mask illumination comprising: an input unit which inputs characteristics of an illumination device; and a processing unit which is configured to change an intensity and shape of an illumination to form an image in an image plane that maximizes the minimum image log slope at user selected fragmentation points.
 32. The apparatus of claim 31, wherein the processing unit is further configured to forcing the intensity at the fragmentation points to be within a predetermined intensity range.
 33. An apparatus for optimizing an mask comprising: an input unit which inputs a desired image pattern; and a processing unit which is configured to change a magnitude and phase of diffraction orders to form an image in the image plane that maximizes the minimum image log slope at user selected fragmentation points while forcing the intensity at the fragmentation points to be within a predetermined intensity range.
 34. The apparatus of claim 33, wherein the processing unit is further configured to obtain optimal mask transmission characteristics by determining horizontal diffraction orders of an optimum mask, wherein the number of horizontal diffraction orders is determined according to the equation: $m = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{x}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}$

where m is the number of horizontal diffraction orders; P_(x) is the pitch of the repetitive cell in the x direction; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.
 35. The apparatus of claim 33, wherein the processing unit is further configured to obtain optimal mask transmission characteristics by determining vertical diffraction orders of an optimum mask, wherein the number of vertical diffraction orders is determined according to the equation $n = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{y}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}$

where n is the number of vertical diffraction orders; P_(y) is the pitch of the repetitive cell in the y direction; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.
 36. The apparatus of claim 33, wherein the step of determining optimum diffraction orders determines optimum diffraction orders in the special frequency domain.
 37. The apparatus of claim 33, wherein the processing unit is further configured to obtain optimal mask by locating areas of maximum transmission and minimum transmission; assigning a primitive area as an area centered on an area of maximum transmission or minimum transmission; and varying edges of each primitive area to match optimal diffraction orders, wherein each primitive areas has a minimum size which is substantially equal to a minimum feature size of the mask.
 38. An apparatus for obtaining an optimum source and an optimum mask comprising: an input unit which accepts user inputs; and a processing unit configured to simultaneously change an intensity and shape of an illumination source and change a magnitude and phase of diffraction orders to form an image in an image plane which maximizes a minimum image log slope at user selected fragmentation points while forcing an intensity at the fragmentation points to be within a predetermined intensity range.
 39. An apparatus for optimizing a placement of transmission and phase shifting features on a mask comprising: an input unit which inputs characteristics of an illumination device; and a processing unit which is configured to obtain optimal mask transmission characteristics based on optimum diffraction orders, locate areas of minimum transmission and maximum transmission, assign primitive areas as areas centered on an area of minimum transmission or an area of maximum transmission, and vary edges of the primitive area to match optimal diffraction orders, wherein the primitive areas have a minimum size which is substantially equal to a minimum feature size of the mask.
 40. The apparatus of claim 39, wherein the mask is a CPL mask.
 41. The apparatus of claim 39, wherein the optimal mask transmission characteristics include horizontal diffraction orders of an optimum mask, and the number of horizontal diffraction orders is determined according to the equation $m = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{x}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}$

where m is the number of horizontal diffraction orders; P_(x) is the pitch of the repetitive cell in the x direction; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source.
 42. The apparatus of claim 39, wherein the optimal mask transmission characteristics include vertical diffraction orders of an optimum mask, wherein the number of vertical diffraction orders is determined according to the equation: $n = {{2\quad {{floor}\quad\left\lbrack \frac{{P_{y}\left( {\sigma_{\max} + 1} \right)}{NA}}{\lambda} \right\rbrack}} + 1}$

where n is the number of vertical diffraction orders; P_(y) is the pitch of the repetitive cell in the y direction; λ is a wavelength of the illumination source; NA is a numerical aperture of the projection optics; and σ_(max) is a radial extent of the distribution of a beam of light from the illumination source. 